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Casimir Forces and Graphene Sheets 
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The Casimir force between two infinitely thin parallel sheets in a setting of N 

such sheets is found. The finite two-dimensional conductivities, which describe the 

dispersive and absorptive properties of each sheet, are taken into account, whereupon 

the theory is applied to interacting graphenes. By exploring similarities with in-plane 

optical spectra for graphite, the conductivity of graphene is modeled as a combination 

of Lorentz type oscillators. We find that the graphene transparency and the existence 

^^ of a universal constant conductivity e^ /{Ah) result in graphene/graphene Casimir 

i-Q interaction at large separations to have the same distance dependence as the one for 

D perfect conductors but with much smaller magnitude. 
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I. INTRODUCTION 

The Casiinir force is a fundamental quantum mechanical relativistic phenomenon which 
originates from the vacuum fluctuations of the electromagnetic field. It couples electrically 
neutral objects with or without permanent electric and/or magnetic moments. In the case 
of two perfectly conducting infinite plates, the Casimir force depends only on the distance 
and two fundamental constants: the Plank's constant and the speed of lighlP^. 

The Casimir force is of particular interest at the nanoscale. It has been shown that such 
a force is responsible for limiting the operation of many nanostructured devices, such as 
nano-electrical-mechanical and micro-electrical-mechanical systems by causing stiction, fric- 
tion, or adhesion^ ^. This has motivated devising experiments to detect the effect in various 
structures. Sophisticated techniques using torsional pendulum or atomic force microscope 
have been used to measure the Casimir force between metallic and dielectric surfaces with 
high accuracy^. The stability of many nanostructured materials, related devices and exper- 
imental settings has also been connected to dispersion forces originating from the Casimir 
effect. Graphitic nanostructures, such as graphenes (a single layer of graphite is graphene), 
carbon nanotubes (cylindrically rolled concentric graphenes), and graphene nanoribbons 
(finite width graphenes) are particular examples^. 

Recently, single layers of graphite have been isolated using micromechanical cleavagd^*^. 
At present, graphene is one of the most interesting and most studied materials and has 
paved the way for future carbon based electronics. Many applications of graphene rely on 
its ability for continuously tuning its charge carrier density and mobility^. This has allowed 
the development of new transistors operating at high frequencies.— Other applications are 
also very promising. Nanomechanical resonators, for example, are especially attractive due 
to their mechanical stability and high resonant frequency^. 

Isolated graphenes also raise the possibility of studying the Casimir force between essen- 
tially two dimensional structures with peculiar dielectric response properties and uncovering 
further insights in the nature of the Casimir interaction. In addition, graphene/graphene 
or graphene/substrate mutual interactions are important components of many experimental 
settings. The Casimir force can be calculated using the Lifshitz theory, which takes into ac- 
count the macroscopic dielectric response of the objects^. This theory was also adapted for 
the case of two graphene sheets using an idealistic description of the dielectric permittivity 



by assuming a Drude type model^. Such an approximation, however, does not take into 
account the electronic properties specific to graphene. Researchers have also considered the 
Casimir interaction between a graphene and a perfect conductor within a quantum field the- 
ory approach by using a more realistic representation of the graphene dielectric response^, 
where low momentum electrons were modeled via the Dirac model. 

The goal of this work is to present a theory for A^ parallel infinitely thin sheets at dif- 
ferent separations and take into account the specific optical properties characterized by the 
conductivity of each separate sheet. The objective is to apply the results to the case of 
parallel graphenes, in order to understand how the distance separation and peculiar charac- 
teristics of their dispersion properties manifest in their mutual Casimir force. Of particular 
interest is the relation of the force to the universal graphene conductivity, which provides a 
different representation of the effect in this peculiar system. Our method utilizes a quantum 
electrodynamical approach based on linear response theory. An essential part is the explicit 
calculation of the dyadic Green's function for this system and the inclusion of the finite 
conductivity of the separate sheets. 

The rest of the paper is organized as follow. In section|TT], the electromagnetic field induced 
fiuctuation stresses between two infinitely thin parallel plates is developed via linear response 



theory. In section III, the dyadic Green's function for two sheets with frequency dependent 



conductivities is found. In section IV, the Casimir force between two parallel sheets within 
N sheets is characterized via recursion relationships for the refiection coefficients. Finally in 
section|V], the theory is applied to find the Casimir force between graphenes. The conclusions 
are described in section IVII 



II. FLUCTUATION FORCES 

In classical electrodynamics, the force is calculated using the Maxwell stress tensor for 
the electromagnetic pressure 

¥ = — Tee - -E^^ + BB - -B^^] , (1) 

where E and B are the electric and magnetic fields, respectively, and 1 is the unit matrix. 
In order to present a formalism to find the Casimir stresses, we consider the radiation-matter 



interaction Hamiltonian within the dipole approximation described by 

5H = - f d^rP{v, t) ■ E(r, t), (2) 

where P is the polarization. 

The quantum mechanical description is performed by simply substituting the classical 
fields with a symmetrized product!^ of quantum mechanical operators (denoted by hats in 
what follows). Then, using linear response theory^, the perturbed electric fields due to an 
external polarization source field become 

dE{r,t)= f dt'^ f d\'?^{r,r',t-t')-P{r',t'), (3) 

^(r,r',t-tO = ^([E(r,t),E(r',t')]), (4) 

where (O) = Tr{pO), p = exp{[FQ — Hol/ksT) is the statistical matrix operator in the 
canonical ensemble, F, is the free energy, ^o is the unperturbed Hamiltonian, ^ is the 
Dyadic Green's function, ka is Boltzmann's constant, and T is the temperature. Note that 
t?(. ., UO ^ t?,r -.,-.) .0. the Kcogeneitv of .pace-«.e. 

Define the structure furrct.ons *?.. and *?bb usmg Fourier transforms as follows 



^ (|E(r, «)E(r', t') + E(r', i')E(r, «)) = f ^*S*EE(r, r', c^)e-"<'-''', 
i {|B(r, i)B(r', i') + B(r', «')B(r, «)) = J ^*^Bfl(i-, r', ^)e-"<'-''>. (5) 



Then using the fiuctuation dissipation theorenpSMl^ one can express the structure function 
f,, iu terms of the .magiuary part of the frequency dependent Green's function at finite 
temperature as 

^EEir, r', u;) = hx ^m^(r, r', u) coth (^) • (6) 

Similarly, the correlations in the magnetic field may be found from Faraday's law by taking 
the curl of the electrrc held. The. are expressed in ter.. of the funct.on V.. by 

VsB(r, r', to) = h^^mV x ^(r, r', cu) x V' coth (^) , (7) 

where xV' is the curl taken on the Green's function with respect to the prime coordinates. 
Therefore, from Eqs.( Tp|[6|[7 ) the temperature dependent quantum electromagnetic stress is 



found via the Green's function^, 



T = — 

An 






i=;£;ir,r ,w 



T 2 = lim / S 

r— !>r' ^ 

Finally, the identity coth(x) = ^' 



27r' 



2 = lim / ^ BB[r,r-,UJ) 

r-!>r' / Zn 



(8) 



n=— 00 x^+n^TT^ 



allows one to obtain 



00 
T 1 = lim ksT y^ G{r,r',iu)n), 

n=— 00 

00 2 

V2 = lim ksT V TT^V X ^{r,r',iun) x V', 



(9) 



where Un = 27inkBT/h are the Matsubara frequencies. 



III. GREEN'S FUNCTION FOR THREE LAYERS 



In order to find the Casimir force from Eq.rtSj) for a specific structure, one needs to 
calculate the Green's function. The Green's function obeys the following equation 



Vx V 



u 



2 
'^(r,r',a;) = 47r^5(r-r'). 



(10) 



which is found from the solution of Maxwell's equations^. 

The system under consideration consists of parallel, infinitely thin sheets located in a 
vacuum. Each sheet is positioned at the boundary between two adjacent layers, and it is 
specified by a two-dimensional, isotropic conductivity 



^ a,,,+i ^ 
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a,,,+i 




(11) 



where j,j + 1 is the boundary between two adjacent layers, cijj+i accounts for the specific 
finite absorptive and dispersive optical properties of each sheet. 
In general, the Green's function may be split in two terms 



^'^- = "di + ^^i 



(12) 
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FIG. 1. The free space regions denoted as 1,2,3 separated by two infinitely thin sheets extending 
in the x — y plane and separated by a distance d in the z-direction. The infinitely thin sheets are 
located at z = and z = d. The conductivities are also denoted. 

where G '-^ is the Green's function in region I with a source in region j, Gg is the free 
Greeu-s fuucfou from a po,ut-l,ke source placed „r layer j without any bouudaries, aud ^'/ 
is the scattering Green's function in region / with a source in region j. The free Green's 
function is later dropped from the calculation of the stress tensor since the Casimir effect 
does not exist in homogeneous space. 

The calculation of the force between planar sheets given in this work is based on a 
procedure using the generalized Fresnel reflection coefficients. The method involves the 
explicit form of the dyadic Green's function for a system of two parallel planar sheets^^'^. 
Thus we fl.t calculate t?(..,., for tire tlrree layer/two .Ireet structure .howrr lu M.0 

Siuce the system has a planar' geometry, we use the general dyadic form of ^ in terms 
of the following orthogonal functions^^l^ 



M(k) = V X [z(; 



and N(k) = |v x M(k), 
k 



(13) 



where = ea;p(zk- r) with k being the wave vector of the electromagnetic excitations. These 
obey the orthogonality relations 



f dVM(k)-N{-k') = 

f dVM{k) ■ M(-k') = {27r)^kl5{k - k') 

f dVNik) ■ N(-k') = {2Trfkl5ik - k') 



(14) 



where k'j_ = k^ + ky and the integration is carried over space. The bulk Green's function is 



found to be of the following form: 

Gl{r, r', u) = — 47r5(r — r')zz + 
{z-z' > 0) 



tu 



lUJ 



dk"^, 



Gr Ur^r',u) = —ATT6(r — r')zz , 
{z-z <0), 



2nklh 

dkl 
27ik^,h 



[M{h)M'{-h) + N(/i)N'(-/i)] 



[M{-h)M'{h) + N{-h)N'{h)] 



J{k±-r±-hz) 



(15) 



(16) 



Mih) = i{ky± - A;,.y)e*^^^-'^+'^^\ M(-/i) = i{ky± - k,y)e' 

M'i-h) = -t{ky^ - k,y)e~'^''^-''^^^''\ M\h) = -z(A;,x - fc,y)e-^('^^-'-'^-'^^'), (17) 

where h = ^JuP' jc? — k\ and r_L = xx + yy . Similarly, N(±/i), N'(±/;,) are defined via 



Eq.(13). 



In this way, the bulk Green's function is expressed in terms of a linear combination of 
transverse electric (M-term) and transverse magnetic (N-term) modes. The scattering part 
of the Green's function is also sought in terms of the orthogonal functions M and N, and it 
is found from the boundary conditions for the tangential electric and magnetic fields across 
each plane. 

The continuity of the electric field across each boundary from FigQ is expressed as 

"^^2(r,r',a;)-^i2(r,r',u;)]=0, (z = 0), 
"^=^2(r, r', uj) - ^22(r, r', a;)! =0, {z = d), 



z X 

Z X 



(18) 



while the tangential components of the magnetic field give rise to surface currents, which 
may be written as 



z X 



z X 



V X ^=^2. / ^ _ y ^ ^ 



12/ 



r,r ,w 



Aniu 



^i,2-^^'(r,r',u;), (^ = 0) 



<y 2,3 



■^32( 



{z = d). (19) 



The surface conductivities account for the finite absorption properties of the infinitely thin 
sheets in the material. The solution of the dyadic Green's function in the different regions 
may be written in terms of the reflection coefficients for the transverse electric (TE) waves 
and the transverse magnetic(TM) waves. For the system considered in Fig. ([I]), these are 

2'nuja2;i/{hc^) _ 27iujai^2/{hc'^) 



Pe 



Pb 



1 + 27iijja2,3/{hc^ 
2Tia2fih/uj 
1 + 2'na2,2,h/u:'' 



Pe 



Pb 



1 + 27ru;CTi,2/(/ic2) 



1 + 27rai^2h/u! 



(20) 
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where (+) superscript defines the top plate and (— ) superscript defines the bottom plate. 

Making the definition Qe,b = 1 ~ Pe bPe b^^*'^'^, the scattering Green's function in region 
two is 

iu"^ f dk'i , 



^ 



22/ / 

'r,r 



C2 7 2TlklW 
-,2ihd 



^M{h)M\h) + ^^^f^ M{h)M'{-h) 

He He 

„+ „— p2ihd ^+ p2ihd 

^PeP^ M(-/i)M'(/i) + ^ m{-h)M!{-h) 

iiE He 

+ ^Nih)N'{h) + ^%^N(/i)N'(-/i) 

„+„— „2j/ia! ^+ p2ihd 

+ PbPb^ N(-/i)N^(/i) + Pj^^{-h)^\-h)]. (21) 



Hb Hb 

Similarly, the scattering Green's function in region one is 



"(r.r' 



' ■ ' c' J 'iTiklfi' 

Z J- (n± J- 9nZnt'lp2iM\ 

M{-h)M'{-h) 



Pe + (Pe + '^PePeV"' 



^E 

+ (PBHpl-2pBPl)e^'-^\ N(_,)N'(-/.)]. (22) 

The other scattering Green's functions are found from the boundary conditions specified in 
Eqs.( p;8|l9| . 



IV. FORCE BETWEEN N PARALLEL SHEETS 

The Casimir force per unit area exerted on each planar sheet Fig. ([I]) is obtained by 
evaluating the zz-component of the difference in the stress between the regions above and 
below that sheet. For example, the force on the bottom one is calculated by taking 

T,= [TfAz)-Tll{z)\_^^^, (23) 

where T^ is the force per unit area on the bottom plate, T^^| is the zz-component of the stress 
in region 2 given a fluctuating source in region 2, and TH is the zz-component of the stress 
in region 1 given a fluctuating source in region 1. The force per unit area on the top plate 
is equal and opposite to that on the bottom one, ie, Tj = — T^. Combining the results in 



previous sections with Eq.(23) one obtains 



T, 



iksT 

27r 



E 



h{iuJn)k±dk_i 



-2ih(iu]„)d 



_p+{iUJn)pE{i^r, 



-2ih{iLj„)d 



_p5(iw„)p^(zu;„ 



(24) 

where h{iun) = i\/{uJn/cy + k\. The expression for Tf, can be used to calculate the tem- 
perature dependent interaction between any two infinitely thin plates in vacuum providing 
that the explicit conductivities are known. One notes that the largest contributiorpS of 



Eq.(24) comes from Und/c ~ 1. Therefore, when hc/{2TxkBTd) ^ 1, the sum above is deter- 
mined by the large n value terms and it can be transformed into an integral with differential 
dujn = 2TTkBT{dn)/h. For T = 0, such a representation is exact, but for T < 300 K, for 
example, d should be less than a micrometer. At distances comparable to the thermal quan- 
tum coherence wavelength At = hc/{kBT) ^ 7 pm, classical thermal fluctuations become 
important and hc/{2TTkBTd) ^ 1 is not valid any more. In this case, the n = and small n 
terms in the sum become important, and the dUn ~ dn transformation cannot be justifled. 
Here, we will assume that we are in a regime where the integral representation is valid. 

The result for the system of three layer/two sheet system can be used to obtain the 
force per unit area in the case of a stack of A^ parallel sheets. Consider the j-th layer in 
Fig. (2). One realizes that the Casimir force results from the inflnite optical reflecting and 
transmitting paths due to the scattering from all the sheets above and below layer j. This is 
described by the effective reflection from below pj^ ^ and from above p^^ b i^ layer j. Then 



combining Eq.(24) with the condition hc/{27TkBTd) ^ 1 and the integral representation 
that follows, the stress on the bottom sheet may be written as 



Tbj 



ih 



27r2 



k\dk\ 



duhiiu) 



-2ih(iui)dj 



PEjii^)pEji^^) 



-1 



-j—2ih{iu))dj 



PBjii^)pBji^^) 



-1 



(25) 



The temperature dependence does not appear explicitly in the force any more. It is ac- 
counted for indirectly through the temperature dependent optical properties of the sheets. 
pf,E,B '^^^ be found via an iterative procedure using a simple recursion relation. Consider 
the three layers denoted as (j — 1), j, (j + 1) by themselves. Due to the inflnite optical paths 



10 

they can be expressed as^^'^ 

PE,j-\,j,j+l = PE,j-l,j + tE,j-l,jpE,j,j+ltE,j,j~ie ' ' + 

iE,j-ljPE,j,j+ltE,j,j-ie ' 'PE,j,j-lPE,j,j+l + • • • , 

oo 

PE,j-l,j,j+l = PE,j-l,j + tE,j-l,jPE,j,j+ltE,j,j-ie'^'^'^' 2^ [pE,j,j-lPE,j,j+ie'^'^'^'] , (26) 

n=0 

where pe,j-i,j is a single sheet reflection coefficient from layer j — 1 to layer j and Iej-ij 
is the coefficient of transmission for a single sheet from layer j — 1 to layer j. Given that 

tE,j-i,j = tE,j,j-i, PE,j-i,j = PE,j,j-i and 1 + PE,j-i,j = tE,j^i,j one obtains 

_ PE,j-l,j + {PE,j,j+l + 2pE,j,j+lPE,j,j-l)e'^'^'^^ , . 

pE,j-i,j,j+i - — — -^Tf^ . [27) 

J- - PE,j,j-iPE,j,j+ie J 

Similarly, one finds the reflection coefficient for the TM modes given that tsj-ij = t^jj-i, 
PB,j-i,j = PB,j,j-i and 1 - PB,j-i,j = tB,j-i,j, 

_ PBJ-IJ + {PB,j,j+l - '2pB,j,j+lPB,j,j-l)e'^' ' , . 

pB,j-l,j,j + l - -^r^ . (.2«j 

J- - PB,j,j~iPB,j,j+ie ^ 
This provides a straight forward method for calculating the reflection coefficients in any 
vacuum layer in a stack of A^ parallel sheets. Suppose j — l,j,j + 1 are a part of the system 
shown in Fig. ([2]). Starting from layer 1, one finds the reflection coefficients between the ffist 
and second layers using Eq.(20). Invoking Eqs.(27p8) recursively by treating the ffist two 



sheets as one, the reflection coefficient in the third layer is found. This is repeated until 
layer j is reached giving the effective reflection from below. A similar procedure is applied 
to find the reflection from all sheets from above layer j, but starting from the top A^ + 1 
layer in Fig. ([2]). 

As an example of the use of the recursion procedure, a four layered system will be con- 
sidered, which corresponds to the bottom four layers of Fig. ([2]). Using Eqs.( 27|28 ), the 
generalized reflection coefficients in layer 2 are expressed as 

+ _ _ PE23 + {PEU + 2p£;23P£;34)e^*'"^^ 

Pe2 — PE23A — ^ 577-7; , 

1 - pE23pE34e^"^'^^ 

+ _ _ PB^^ + (PB34 - 2piJ23PiJ34)e^*^'^3 

Pb2 — PB234 — ^^TX , 

1 - pB23pB34e^""^^ 

Pe2 = PE2U Pb2 = PB2l (29) 



Note that Pe b2 ^^^ actually p^ ^ from Eq.(20), since the system below layer 2 is the same 



as the one from Fig.(l). The reflection coefficients can be substituted in Eq.(25) yielding 



the Casimir force per unit area on the bottom plate in the four layer /three plate system. 



11 



N+1 



L 



(J 

N,N+1 ■ 








7., 








■ N 




• 


N 






^j,j*1 ■ 


■ 


j+1 




7 






J 


^j-1,j ■ 




i ^i 


■^M 


■ 


j-1 




^2,3 . 


■ 


3 




, 7 








^2 






2 




■ ^1 


^1,2 ■ 




1 





FIG. 2. A^ infinitely thin sheets located in free space and separated by distances dj. The sheets 
extend in the x — y plane. Their conductivities and positions along the z-axis are also shown. 

V. CASIMIR INTERACTION BETWEEN GRAPHENES 



Before the graphene planes are considered, it is useful to calculate the interaction for the 
limit of infinitely conducting planar sheets. In that case, the conductivity becomes cr — > cxo 



yielding p^- — t- — 1 and p^- — t- 1. Thus using Eq.(25) we recover the well known resultPfor 



the magnitude of the attraction between two parallel perfectly conducting plates separated 
by a distance d 

\To\ = TTTTT^I- (30) 



240^4 



A. Universal conductivity 



Further, we apply the results for infinitely thin sheets obtained in section IV to calcu- 
late the Casimir force between two graphenes. Researchers in the past have considered 
graphene as an infinitely thin sheet, and have shown that this is a reasonable approxima- 
tion for distances greater than a few times the interlayer graphite separation^SHSEZI ^]iq 
emphasis now is to specify the graphene conductivity. It has been predictecP^'^ and found 
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experimentally^, that over a relatively wide range of photon energies (up to 3 eV), the 
graphene conductivity is approximately constant given by the value of ctq = e'^/{4h). This 
peculiar effect is closely related to the energy band structure of graphene, as shown in the 
Appendix. 



Given a constant conductivity ctq for two parallel sheets separated by a distance d, Eq.(25) 
may be written in the following form 

3hc ^fl 



167r2# 



n" 



[F{ao,n) + G{ao,n)] 



F{ao,n) 
G{ao,n) 



n=l 

c ^, 27rcro/c ^ . .\ 
/3( , , ^ , ,2n + l,-l). 



27iaQ 1 + 27ro"o/c 
27rcro / 1 



2n- 1 



where 



(3{x, a, h) 
is the incomplete beta functiorpi'. 



1 



e~Hi 



27rcro/c 
1 + 27rcro/c 

- tf-Ut 



(2n-l)' 



(31) 



(32) 



For graphene, however, 27rcro/c <^ 1. Then Eq.(31) is approximated by taking the first 
term in the sum, which reduces to 



IT I ^ 



"ihc 

87T^d^ 



1 - 



27ro"o 



ln(l + 2nao/c) 



Further expansion of the In-function and inserting o"o gives l^gl ^ lei^^T^ 



3e^ 
327rd2 



(33) 



. Thus 



the leading term in the force does not depend explicitly on the Plank's constant and speed 
of light any more. This is a remarkable result originating from the particular value of the 



graphene conductivity. We note that the approximate result from Eq. (33) is fairly accurate 



since it differs by less than 2% from the numerical integration of the exact result (Eq.(25)) 



It is also interesting to see the similarities and differences between Tg and Tq. In particu- 
lar, the distance dependence of the Casimir force is the same as the one between two perfect 
conductors. However, comparing their values gives Tg/To ^ 0.00538. Thus the graphene 
Casimir interaction is much smaller in magnitude than the interaction between perfect con- 
ductors. This is directly related to the transparency of the graphene system, reflected in its 
small constant conductivity value. The interaction between a perfectly conducting plate and 



graphene can also be calculated via Eq.(25). In this case the force is much larger as com- 



pared to the one between two graphenes. Given the constant ctq for graphene and a — )■ 00 
for a perfect metal, one finds T/Tq ^ 0.025. 
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The general formula from Eq. (25) allows one to consider the Casimir interaction in 



a system with more than two parallel graphenes. We can estimate the force between two 



graphenes in a three graphene setting. Using Eq. (25), the force is obtained as 



'^' = 2;^ 16^/ ^/o ^'^^E^~'"[(/^^2lPi.234r + (pB2lPi.234r], (34) 



where the reflection coefficients were defined in Eq.(29). Since 27r ctq/c ^ 1, the dominant 
contribution comes from the n = 1 term: 



?)hc 



2nao 2 /27rcro V 



do 



^2 + ds 



(35) 



IGTT^t/l \ c 3 \ c / 

Thus the interaction is determined mainly by the two adjacent graphenes in the three 
graphene system, and it is affected little (only to order ctq) by the presence of the third 
one again due to their transparency. 

B. Other models for the conductivity 

The low energy graphene band structure has been very successful in explaining experi- 
mentally observed properties at various temperatures^^'^. Since the Casimir interaction at 
larger separations is determined by that regime (corresponding to optical excitations less 
than 3 eV), one concludes that its qualitative features cannot be an exception. Nevertheless, 
as the graphenes are brought closer, the presence of the higher energy bands besides the 
ones closest to the Fermi level also needs to be considered. 

Direct measurements of the graphene conductivity have been done at photonic energies 
below 3 eV so far^. Experiments appropriate for higher regimes have not been currently 
reported. At the same time, investigating a theoretically is one of the issues at the fore- 
front of graphene science. Recent ab initio calculations indicate that the in-plane optical 
properties of graphite and graphene are very similar over a wide range of frequencieP^ES. 
Experimentally it was also reported that the optical conductivity of graphite per graphene 
sheet is very close to the universal ctq value of an isolated graphene^ for low optical frequen- 
cies. Thus a viable approach for further investigating the Casimir force between graphenes 
is to use the in-plane optical data of graphite and transpose it to graphene. 

Results from ab initio calculations for graphite have been mapped to a series of Lorentz 
oscillators with a Drude term, whose parameters fit previous graphite measurementsSS be- 
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FIG. 3. (Color online) a) Re and Im parts of the dynamical conductivity for in-plane graphite. 
The universal graphene conductivity do is also shown. The insert displays the infrared regime. 
b)The Casimir force between graphenes as a function of distance with constant and dynamical 
conductivities. The force is normalized to the one for perfect conductors. 



tween 0.1 eV and 40 eV. A plot of the in-plane conductivity as a function of photon energies 



is shown in Fig. (3(a) ). We note that in the infrared spectrum a of the two systems are dif- 



ferent - insert of Fig. (3(a)). For graphite, the conductivity exhibits a Drude-like behavior 
from intraband transitions, while the conductivity for graphene stays constant. This has 
also been observed experimentally^ and it can be explained in terms of the electronic struc- 
tures of the two systems. The onset of the Drude-like term in graphite is related to the 
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splitting of the energy bands and their becoming shghtly parabohc due to the interlayer in- 
teraction, while the constant o"o in graphene originates from cancelations occurring between 
the intraband and interband transitions due to the linear in k energy bands as shown in the 
Appendix. Thus for the calculations here, we modify the fitted model^^ for photon energies 
in the infrared region (below 0.05 eV) by requiring a = ctq as displayed in the insert of 



Fig.(3(a)). 



For larger photonic energies, the in-plane graphite conductivity is mainly determined 
by the single graphene properties. The Lorentz oscillator model shows that Re{a) stays 
relatively constant in the low optical regime (up to 3 eV), which means that the univer- 
sal graphene conductivity has not been affected significantly by the presence of the other 
graphene layers. This is in agreement with previous experimental findings^. Also, the two 
peaks in Re{a) that appear around 5 eV and 15 eV are related to tt — vr* and a — a* elec- 
tron transitions for an isolated graphene^^l^, respectively. A sizable imaginary part of the 
conductivity also appears after 3 eV. 



Using this model, the force between two parallel graphenes is found via Eq.(25). We 
show a plot of the normalized to perfect conductors force per unit area as a function of 
distance in Fig.(|3(b)) with a for graphite with the modification of cr = ctq for energies less 



than 0.05 eV and for a = cxq over the entire range. The plot shows that at longer distances 
the force approaches the one given with a constant conductivity, but at shorter distances 
higher photon frequency modes contribute to an increasingly larger force. 

VI. CONCLUSION 

We have studied the Casimir force between parallel infinitely thin sheets in free space. 
The particular absorption optical properties are taken into account via the conductivities 
of each sheet. The derived expressions rely on generalized Fresnel reflection coefficients ob- 
tained with an iterative procedure using the Dyadic Green's function. This is especially 
convenient since it is applicable to a system of A^ sheets. The theory is applied to the 
case of graphene/graphene Casimir interaction in order to study how this fundamental ef- 
fect depends on the graphene optical response and the distance separation. The graphene 
conductivity is described with a model based on the low energy band structure first. In 
this case, we find that the Casimir force obeys the same distance dependence as the force 
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between two perfect conductors, but it is much smaller in magnitude due to the graphene 
transparency. These results are directly related to the existence of a constant graphene 
conductivity cq = e^ /{Ah) over the optical range of photon energies. This universal value 
translates into a mutual Casimir force that depends only on the electron charge and the 
distance. Because of its transparency, the interaction between two graphenes is not affected 
significantly when more graphenes are present. 

The graphene conductivity is also calculated using a model based on ah initio calculations 
and appropriate for in-plane graphite optical data due similarities between the two systems. 
For graphite, a however, has to be modified in the infrared photon energy region in or- 
der to reflect experimental and theoretical results for the existence of a constant graphene 
conductivity. This is important for the interaction in the limit of large separations. 

Finally, we comment that the graphene conductivity might be influenced by other factors 
such as electron-phonon interaction, electron correlation effects and the presence of exciton. 
Efforts to understand and quantify these effects theoretically have just begurP^. Experimen- 
tal research is also needed to validate such studies. Thus it would be interesting to explore in 
the future other models for the conductivity of graphene in relation to their mutual Casimir 
interaction. 
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Appendix A: Conductivity Model 

The conductivity of graphene can be modeled using the low energy electron excitations 
which obey a linear momentum energy dispersion relation e = ^vpk, where vp ~ c/300, 
and k is the magnitude of the two-dimensional wave vector^. Within the Kubo formalism, 
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FIG. 4. (Color online) Graphene conductivity a{iuj) in units of (Tq = e^/(4/i) vs. frequency is given 
at two temperatures, T = 300 K and T = 30 K. The scattering rate is F = 0.1 eV. 

the conductivity is expressecF^'^ using 



O'intrai^, T) 



Cinteri^, T) 



'nh?'{uj + iF) Jq 
ie^fw + zF) 



erfe 



C^(^5 F) = aintrai.^, T) + ainteri.^, T), 

^e^ r . fdf{e) df{-e) 

de de 

_ /(-^) - /(^) 

7rn2 ./n ""^u + iVy-A^e/hy 



de- 



(Al) 



where /(e) = l/[exp(e//i;BT) + l)] is the Fermi-Dirac distribution function and F is a damping 
parameter which accounts for physical processes contributing to broadening of the optical 
spectrum, aintrai'^-,'^) is the intraband contribution to the conductivity which is found to 
be 



mtray , > -nh^UJ + iV) ' 

while the intraband contribution to the conductivity is calculated as 



(A2) 



(^intray^-i T) 



SnkBT 



tanh(x) 
ax 



x^ 



ikgT 



2- 



(A3) 



In the small temperature limit, A2 and A3 result in a being a constant - ctq = e^/(4/i). 
Our sample plot of a{iu, r)/o"o vs. u, however, shows that even at higher temperatures, the 
conductivity does not differ much from the universal constant value. Our calculations also 
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show that the temperature entering through the conductivity has httle effect on the Casimir 
force. Indeed, the force at 300 K differs by less than 1% as compared to the one for OK. 
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